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HEAT TRAESEEH I2T TEE TuRBUIjEET BOUEDARY XAYSE 
0? A COMPRESSIBLE G-AS AT HIGH SPEEDS 1 
By E. Erankl 



SUMMARY 



The Reynolds law of heat transfer from a wall to a 
turbulent stream is extended to the case of flow of a 
compressible gas at high speeds. The analysis is "based 
on the modern theory of the turbulent boundary layer with 
laminar sublayer. The investigation is carried out for 
the case of a plate situated in a parallel stream. The 
results are obtained independently of the velocity dis- 
tribution in the turbulent boundary layer. The heat 
transfer q. in unit time is a function of the frictional 

coefficient \1» = — — , the free stream temperature f , 

pus 

the temperature at the plate T* , the free— stream 

velocity u and the density p , the velocity at the 

boundary of the laminar sublayer and the basic turbulent 
layer u, , and the Prandtl number for the given gas 




The problem of heat transfer at high velocities has 
at the present time become a very pressing one both in 
connection with heat technology and high— speed aviation. 
Recently this problem has been encountered in the computa- 
tion of the surface of a wing radiator. The existing 
papers on heat transfer f r os the wall to the gas in part 
refer to small speeds (Reynolds, Prandtl, and others) 
while with regard to high— speed flow some ( Bus em an , 
Erankl) consider only the laminar regime, others (Crocco) 
neglect the laminar sublayer in the turbulent layer, 
while still others (G-uc hman, Shirokov) start out from 

iReport Eo. 240, of the Central Aer o— Eydr odynamical Insti- 
tute, Moscow, 1935. 
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relatively rough hydraulic assumptions. In the present 
paper the problem is treated from the point of view of 
the aodern turbulence theory. 

We shall consider the case of plane parallel flow 
and in particular the boundary layer of a plate in a 
stream parallel to it. We shall first review some basic 
known r esult s . 

The instantaneous velocity components u, v and 
the values of the state variables p, p, T are broken up 
into their mean values u, v, ... and fluctuations about 
the mean values u 1 , v* , ... We assume* that the turbulent 
fluctuations of the density p' are small by comparison 
with the mean density p. We then have the following 
equations for the turbulent stresses T xx , T xy» T yy and 
the heat transfer per unit time through unit area q. :o o_yi 
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iPhe viscosity forces are here neglected by comparison with 
the turbulent friction. 

We now proceed to the case of the plate. In this 
case p has a constant value. We choose the edge of the 
plate for the x— azis. We may then neglect the deriva- 
tives with respect to y, and v by coaparisor. with u. 
We assume, finally, that the average energy of the flr.c— 
tuating turbulent velocities has a constant value. Equa- 
tions (o) and (5) then "become 

^ br *Z (3.) 



dt dy 



dt V 



'°p i + ¥)--& ( * t ^- ,2 £ (5,) 



3For the determination of T X y and q.^ we start 

from the following generally accepted assumptions. 2he 
deviating of a "body o-f fluid arriving at the point (zy) 
from the point with ordinate y + l y , is given by 

» 1 du , „ s 



Similarly we obtain 



7 dy 



J oy 

Substituting these expressions in (l) and (2) we obtain 



where A is the coefficient of turbulent apparent vis- 
cosity, A the coefficient of turbulent heat conductivity. 
From the equations 



A = pr«l y 
A = c pP v' l y 
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there follows directly the well-known equation of Prandtl 

A = c,- ; A ( 10) 

Substituting equations ( 8 ) and (9) in equations ( 3 1 ) and 
(5 r ) and dropping the averaging signs, since in what 
follows we shall deal only with mean values, we obtain 

du 5 f ou "\ . . 

d / us x a r a / uA] , . 

' dt V ? 2 J by I by V * 2 J J 



Thus the velocity and total energy Jc T + — 

- it 

satisfy the sane linear differential equation. All bound- 
ary conditions may therefore be satisfied by setting 

Jc. 0 T + ~ = an + b (13) 

where a and b are independent of y. This relation 
was first derived by Crocco, 

Yfe now proceed with our problem, requiring from the 
above relations only equations (10) and (13). We shall 
use the following notation: 



e = J c v T + ~- total energ 



u, T , e, . . .the values of u , T , e,...in the free flow 
vl x , T x , e , the values at the boundary of the laminar sub- 
layer and turbulent layer , T* , p* , the values of T, 
p, p. at the wall. 

Equation ( 13 ) may then be rewritten in -the form 

e — e = 



e — e 



1 T7 _ - l 



1 (u - u.) (14) 



In the laminar sublayer the tangential stress T,. y 
and the heat transfer q., r are determined by the formulas 
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xy 
1v 



du 

* SI 
*7 



where jj, is the coefficient of viscosity and \ the 
coefficient of heat conductivity. Prom the physical 
significance of t xy and q. y it follows that they 

should he continuous over the entire region of flow. 
(The schematic representation of completely laminar and 
completely turbulent sublayers leads therefore to a dis- 

o" u ^ 

continuity of ^ and — in passing from the laminar 

to the turhulent layer since \t is considerably less 
than A and X less than A). 



As a first approximation we assume a 
hut ion of velocity and temperature in the 
layer. Ve then have 



linear distri- 
laminar sub— 



u x = 



T 



(15) 



T. 



(16) 



where t and q are the values of T 



and 



xy ailu - 1y 

directly at the wall and h is the thickness of the 
laminar sublayer. Prom the assumption of the continuity 
of T^, r and q y it then follows that 



' xy 



5u 
•i "oy 



= T 



(17) 



. St 

* dy 



= a 



(18) 



Prom equations ( 10) , (17), and ( 18 ) we obtain the 
relation 
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q, du 



(19) 



We now differentiate equation (14) with respect to 
y for y = h and obtain 



he 



e — du 
u-u, dy 



(20) 



Substituting the value ^— ! from (13) we obtain 

6y \i 



Ja e — e, 

u + — - = — 

1 T U — U 



(21)- 



Finally substituting the value of H ± , from (15) 
and (16) 



X l ** X* T U l 



(22) 



there is obtained the following relation between q, t, 
u, T, T* and u 1 : 



Jc 



Jq 



P 



<1 J*, 



\ 



U"= u 7 * 
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u. + 
1 T 



u — u. 



Solving this equation for q we obtain 



c p (T - T ; 



) + (^_z_iil: 



= T 



u + u 1 (Pr - l) 



(23) 



r - ' i* + 



(u - u,) : 



2J c. 
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where Pr = — jj— (= 0.8 for diatomic gases). 1 

The heat transfer through unit area per unit time 
can therefore be found when the temperature in the free 
flow and at the wall, the velocity of the flow, the 
coefficient of friction, and the density at the wall are 
known. The velocity u z on the right side of the equa- 
tion can he determined in the following manner. In an 
incompressible fluid according to the hypothesis of 
Earman we have 



(25) 



where s is an absolute constant. According to the 
tests of Kikuradse it has the value 

s = 11.6 (26) 
In a compressible gas we have the analogous relation 

where s may in general depend also on the nondimens ional 
parameters — and 



P X*T* JTpl 

S ince , however, s enters only into terms that are 
comparatively small at not too high (in particular, sub- 
sonic) velocities we may as a first approximation use the 
value s = 11.6. Prom (25') we have 



A* 




~2~ (27) 



so that equation (29) may be written also in the form 

1 In all our formulas Cp denotes the heat capacity of a 
unit mass at constant pressure. If unit weight instead of 
unit mass is used then Cp is replaced by gCp. 
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q. = 



- T, 



"5^ 

2y 



_U_3_ 

2JC-P 



( 25 



1 + (Pr - l)s'V ^* 



C0HEECTI01T 



These papers had already gone to press when it was 
pointed out to the author "by II. Shir okov that according 
to the unpublished results obtained "by him it is not 
correct to neglect the dissipation of energy in the 
laminar sublayer. As a result the curve of temperature 
distribution in the laminar sublayer cannot be considered 
a straight line. Instead is obtained the parabola 



= T* 



7 s 



Equation (23) above is therefore replaced by the 
equation 

. - s u 3 + (Pr - l)u s 

c(T* — T ) L - 

P * 2J 
a _= t _ 

u + (Pr - l)u 1 

and. similarly equation (12) in the second article is 
replaced by the equation 

i = l + JL r cp + s (Pr - l)] - v[cp s - (Pr - 1) s s l . 

CT Pr 2i-c 

Correspondingly there is a change in the remaining 
conclusions of our articles: in particular, there is a 
quantitative change in the dependence of the friction 
coefficient on Ba. We may also note that A. G-uchman 
and his coworkers (Leningrad) found experimentally a 
considerable lowering in the friction coefficient in 
pipes with increasing 3a. 
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r£ICTI0H II? THE TTJBBTJLEST BOTJSDASY LATE2. 
OJ A COMPHESS IBLE GAS AT HIGH SPEEDS 
By 1", Prankl and V. Volshel 



SUMMARY : 

The Prandtl— Xarman theory of turbulent friction, the 
so-called "Mis chungsweg" theory, is extended to cover the 
case of high— speed motion of a compressible gas with heat 
transfer from the wall. A correction of the Sarin an fric- 
tion law is derired that depends on the Bairstow nunber 
(the compressibility effect) and on the ratio of the abso- 
lute temperatures in the free stream and at the wall. The 
analysis is based on the heat transfer law developed in a 
preceding article by 5". Erankl. 

The airfoil drag is made up of two parts, namely, the 
drag of a plate of the same chord and the "form drag," that 
is, the excess drag above that of the plate. The effect 
of the compressibility on these two drag components depends 
on a number of factors. The form drag coefficient varies, 
mainly with change in the velocity distribution about the 
airfoil while the plate drag coefficient depends only on 
the change in the friction coefficient. In the present 
paper we shall investigate the effect of the compressibil- 
ity on the frictional coefficient and also the effect of 
heating the wall, having in mind the case of the wing 
radiator. We restrict ourselves to the case of the turbu- 
lent boundary layer of a plate situated in a flow of con- 
stant velocity at zero angle of attack. Eor this purpose 
we generalize the new theory of the turbulent boundary 
layer, the so— called mixing path ( His chungsweg ) theory. 

Before writing down the fundamental differential 
equations we pass to nondimens ional magnitudes by intro- 
ducing measures for the lengths, velocities, and gas state 
variables. The nondimensional velocities and state vari- 
ables as functions of nondimensional coordinates will 
then depend on parameters which are nondimensional combi- 
nations of the measures, the-' boundary conditions at the 
plate, and the physical constants of the gas. As boundary 
data we may take the frictional stress T the heat trans- 
fer per unit time q. and the' state variables of the gas 
p, p # , apd T„ at the plate* The pressure p, as 
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usual in "the "boundary layer theory, is assumed constant. 
We now introduce p„« and as measures for 

As a measure of the velocity we shall take the 
velocity" 



p and T 
" f r i ct ion 



v * = 



(1) 



and as measure of the length 



P* v * V P*t 



(2) 



vfhe're p. ^ is the value of the viscosity coefficient at 
the plate. 3y introducing the above measures the param- 
eters of the solution may now he only nondim ens i onal com- 
binations of the "boundary data and the physical constants 



di + 



Dgether with Ear man ,. ........„.-.., 

" J --J-J-- near the wall 



.. _ we assume that 

tions at points very ' ~ 



nat only ttie end c 
produce an effect . 



the end con- 



In the absence of compressibility and ho at transfer 
such combinations do not exist, as a result of which fact 
there is obtained as is known a universal velocity dis- 
tribution not dependent on any parameters. We shall 
therefore in what follows denote these measures as "uni- 
versal." The only such independent parameters are 



T 



(3) 



* = 



**** V P*T 



(4) 



and the magnitudes depending only on the physical constants 



P.r = 



• c -oM- 



(5) 



(6) 



Jt 
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For diatomic gases we have k = 1.4 and the Prandtl 
number Pr = 0.8. (In the above c p , c v are the specific 

heats at constant pressure ana volume, respectively, \ 
the heat conduction coefficient of the gas, X + the value 
of X at the wall.) 

We shall now investigate the turbulent boundary layer 
of the plate. ¥e shall take the edge of the plate for the 
x— axis; the variable velocity in the boundary layer will 
be denoted by u; the other notation is the same as for 
the preceding article "Heat Transfer". 

We now assume that for the compressible g*as the 
fundamental differential equation of the Zaraan turbulence 
theory remains valid: 



K 2 p = T (7) 

W 3 > 



where K. is an absolute constant equal to 0.4 and t is 
the frictional stress at a distance y from the plate. 
When the coordinate x (that is, the distance from the 
forward edge of the plate) is sufficiently large, however, 
we nay neglect the forces of inertia. 3 e cause of this and 
from the constancy of the pressure p it follows that T 
does not depend on y. 

We now proceed to the "universal" measures and in- 
troduce the nondimensional magnitudes; 



tj = 



o- = 



JL 
p* 



u 

h 

T 



(8) 

(9) 
(10) 



Equation (7) then assumes the form 
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Vor, 2/ 

7 = 



In order to obtain the differential equation for cp 
it remains only to express a in terms of' cp.. This is 
possible with the aid of equations (14) and (2l) of the 
preceding article. Prom these there is obtained an equa- 
tion analogous to equation (23) of the preceding article: 



q = T 



/ N - (u - u. ) 

c I T - I ;+ i- 

P \ J 2Jc p 

u + u (Pr - l) 



which after substituting the universal measures is re- 
duced to the form 

_ 1 v _ -, 

i = 1 + — 9 + s (Pr - 1) I - 7 (9 -s) 3 (12) 

cr Pr «_ J 2k ' 

where s = u 1 /v # according to the above considerations 
can depend only on Y, *, k, and Pr . 1 

Substituting expression (12) in (ll) we obtain the 
differential equation for cp which is easily solved in 
terms of elementary functions. 

To determine the boundary conditions we start from 
the hypothesis of linear distribution of u in the 
laminar sublayer. This gives 

T 

u = — y 

or 

<P = ^ (13) 



and hence 



cp 1 = cp ( s ) = s 



1 See Correction on p. 8. 
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As regards the value of ccp/dTi in the turbulent 
layer for y = s (that is at the boundary of the laminar 
sublayer) we have 



&9 

5n 



= f 



(14) 



n = s 



where f, like s, can depend only on Y, «, k and Pi 
In an incompressible fluid in the absence of heat trans- 
fer f has the constant value 



f = 0.289 



(15) 



Our problem now is to find the coefficient of fric- 
tion \|/* (see preceding article) as a function of the 
Bairstow number, the ratio of the temperature in the free 
stream to that at the wall and the Seynolds number at a 
given point of the plate, that is. 



I. 3a*. 

5 J 



(16) 



where 



Ba* 



u 
a* 



(17) 



a* = «/(ir - 1) Jc p T* 



(= a* velocity of sound at the wall) 

Instead of T^/T it is convenient in what follows 
to introduce the parameter u> determined by the equation 



f 



k - 1 



Ba^u) + 1 



(18) 



In what follows we shall assume that w is of the 
order of magnitude of unity. The physical significance 
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of this assumption lies in the fact that T* - I is of 
the same order of magnitude as the temperature differ- 
ence that occurs with friction and without heat transfer. 
Jor wing radiators at large velocities this assumption 
may "be considered as satisfied. 

It is now necessary to determine %'* in the form 
\|/* = 3a* , u>) 

lor this purpose it is necessary first of all to ex- 
press the parameters 7 and v in the velocity distri- 
bution, law in terns of \|/ * , 3a^ and . uu . 

lor y we obtain immediately 

Y = j- 3a* 2 (19) 



and for •& from equation ( 38 ) of the preceding article 



2 



* - 



k - 1. 



3a, 



1 + (?r - l)s 



UD - ■ 1 



i (20) 



¥e now substitute ( 19 ) and (20) in the solution of 
equations (ll) and (12) for. the boundary conditions (13) 
and (14). If this solution is expanded in a Caylor series 
in powers of 3a,,, 3 and only the zero and first power 
terns retained there is obtained 



,K(cp - s) 



T| = S + 



K f 



4K 2 f I 



X) I I (m - cp) 2 + n 



— r r 



(21) 



where m, n and r depend only on w , Pr , 'I'*, k . and s . 
In the case 3a* = 0 we obtain the well— known logarithmic 
velocity distribution. 
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Por the solution of this problem it is still re- 
quired that the velocity distribution in the well— known 
integral condition of Xarman he substituted: 



4 
dx 



Jpu ( u — u ) dy = 



= — T 



or in absolute measure 



d.i - - ~- 0. i f aw ( qj - a? ) ±r\ 
v# i i 



where 



! o 



J 



(22) 



Carrying out the integration of equation (22) on 
the assumption that' f and s maintain the constant 
values s = 11.6 and f - 0.289 and setting 



we obtain 



I = P 1 (z) + e s ? £ (z) - 3 a «a^S 0 (u>) + S^kOz + S a ( n ) ln(z- c ) 



+ C In z + 



;( z — c) : 



+ S 6 (c) z + S 7 (ta) z 2 + 



C 3 + S 8 {m)z + S s (o)gg 

Z( B - c) 3 



+ C4B3.Cz) + S 10 (to) Si(z - c)> 



( 23) 



where 



P a (z) = 9542 - 3134z + 245: 



P„ ( z) = 1.31 ( z^ - 4z + 6) 
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S 0 (u)\~ 86217 - 2675 (<u - l) + 168(u)-l) 2 



s x ( 


a) 


= 129.7 (a - l) - 33.0..(u)-l) 2 
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(0 — 1 ) 3 


] 


c = 


0 


.928 










13 433 






c 

p. 




136231 










- 58. 416 






c 4 




18 4.94 








O 


z 

f B 7 ' Cl 7, 
1 

J Z 

— oc 






We 


th 


•as have an equation in the form % 







and our fundamental prohlem is solved. It regains to not 

only that the dependence of the Heynoldn ranker 3.eg - 

p TES/u^ on \|; # , Ba # , and u) can likewise he readily fen. 

Zeiss , 
C OHa it t ee 



Translation by S. 
National Advisory 
for Aeronautics. 



